Chlorine Adsorption on Graphene: Chlorographene 
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We perform first-principles structure optimization, phonon frequency and finite temperature 
molecular dynamics calculations based on density functional theory to study the interaction of chlo- 
rine atoms with graphene predicting the existence of possible chlorinated graphene derivatives. The 
bonding of a single chlorine atom is ionic through the transfer of charge from graphene to chlorine 
adatom and induces negligible local distortion in the underlying planar graphene. Different from 
hydrogen and fluorine adatoms, the migration of a single chlorine adatom on the surface of perfect 
graphene takes place almost without barrier. However, the decoration of one surface of graphene 
with CI adatoms leading to various conformations cannot sustain due to strong Cl-Cl interaction 
resulting in the desorption through the formation of CI2 molecules. On the contrary, the fully chlo- 
rinated graphene, chlorographene CC1, where single chlorine atoms are bonded alternatingly to each 
carbon atom from different sides of graphene with sp 3 -type covalent bonds, is buckled. We found 
that this structure is stable and is a direct band gap semiconductor, whose band gap can be tuned 
by applied uniform strain. Calculated phonon dispersion relation and four Raman-active modes of 
chlorographene are discussed. 

PACS numbers: 73.22.Pr, 63.22.Rc, 61.48.Gh, 71.15.Mb 
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I. INTRODUCTION 

The synthesis of graphenei 2 - has led to intense re- 
search activity in the field of graphene-based nanoscale 
devices. Although graphene is one of the most mechan- 
ically strong material having a wide range of extraordi- 
nary properties^— practical device applications are lim- 
ited by its metallic behavior and sensitivity to surface 
adsorbates. 

Efforts to synthesize chemically modified graphene 
composites with tailored electronic, optical, and chemical 
properties have presented new directions in graphene re- 
search. In particular, band gap engineering of graphene 
through chemical modification, such as oxygenation,^— 
and hydrogenationi 2 - - — is appealing for electronic appli- 
cations, since the scalable fabrication of graphene-based 
devices without disturbing the strong honeycomb lat- 
tice has become possible. However, due to the complex 
atomic structure of grapheneoxides 7 - (GOs) and thermal 
instabilities of hydrogenated graphenes (CHs) even at low 
temperatures,— ~— search for the novel graphene-based 
materials is still continuing. 

Over the past three years experimental^— and 
theoretica l 22 ' 23 studies have demonstrated that chemical 
conversion of graphene to fluorographene (CF) is possi- 
ble. In addition to early studies on the atomic composi- 
tion and band structure of fluorocarbon materials ) 24 ' 25 it 
was reported that the monolayer CF has quite different 
vibrational spectra and Raman characteristics as com- 
pared to hydrogenated graphene analogues.— We also 
investigated the electronic and elastic properties of pos- 
sible fluorinated graphene materials and attempted to 
clarify the discrepancy between theoretical and experi- 
mental results^ Easy synthesis, high-quality insulating 



behavior and extraordinary mechanical strength of fluo- 
rographene (CF) have inspired intense research on other 
halogen decorated graphene derivatives. 

In addition to three known derivatives of graphene: 
GO, CH and CF, the successful synthesis of chlorinated 
graphene was also achieved very recently^ It is experi- 
mentally demonstrated that nondestructive and pattern- 
able conversion of graphene is possible by using var- 
ious photochemical chlorination techniques i 2 ^— While 
the research on chlorine-graphene interaction is rapidly 
growing,^ - — comprehensive research on the stability of 
various chlorinated graphene structures and their result- 
ing properties are sparse. 

In this paper we present a detailed analysis of the in- 
teraction between chlorine atom and graphene leading 
to the chlorination of graphene. Although the possibil- 
ity of covering graphene surfaces with chlorine atoms has 
been reported, analyses of structural stability, electronic 
and magnetic properties as a function of CI coverage are 
lacking. Our main motivation is to reveal which con- 
formations of chlorinated graphene are stable and how 
these conformations modify the properties of graphene. 
To this end we investigated the chlorination of graphene 
starting from single CI adsorption to full coverage lead- 
ing to chlorographene, namely CC1. At low coverage 
with diminishing Cl-Cl coupling, the binding of CI to 
graphene is significant, but adsorbed CI atoms migrate 
on the surface of graphene almost without an energy 
barrier. We found that the nonbonding chair conforma- 
tion of chlorographen o 30 ' 31 (CC1) consisting of a planar 
graphene sandwiched between widely spaced two planar 
CI layers is unstable. On the other hand, the covalently 
bonded chair conformation of the chlorographene (CC1) 
is found to be stable at T=0 K and possibly at room 
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temperature. This latter conformation consists of buck- 
led graphene sandwiched between two planar CI layers 
and is a nonmagnetic semiconductor with 1.2 eV direct 
band gap. Our results reconcile the discrepancy between 
the experimental study2£ obtaining semiconducting prop- 
erties upon the chlorination of graphene and theoretical 
studies predicting metallic statei 30 ' 31 



II. COMPUTATIONAL METHODOLOGY 

To investigate mechanical, electronic magnetic prop- 
erties of chlorinated graphene we carried out first- 
principles density functional theory (DFT) calculations 
within the local density approximation (LDA)2£ using 
projector augmented wave (PAW) potentials^. All re- 
sults discussed in the text are obtained using LDA. To 
compare with the LDA results of specific systems, we 
also performed calculations using Generalized Gradient 
Approximatior>22,(GGA) together with van der Waals 
(vdW) correction, (GGA+vdW) 4£ In our earlier tests 
LDA yielded interlayer spacings of layered materials and 
other structural parameters in agreement with exper- 
imental data, as well as with those obtained by us- 
ing GGA+vdW. In the vdW corrections of the later 
method, DFT description is restricted to shorter cor- 
relation length scales, but for the medium and large 
inter-atomic distances the damped C§r~® term is used. 
The systems, whose numerical values are also obtained 
by using GGA+vdW are indicated. Numerical calcula- 
tions are performed using VASP4i Kinetic energy cutoff, 
E C ut = ?i 2 |k + G| 2 /2m, for plane-wave basis set is taken 
as 500 eV. A vacuum spacing of at least 15 A is placed 
between adjacent layers to hinder the layer-layer interac- 
tions. The convergence criterion of self consistent calcu- 
lations for ionic relaxations is 10~ 5 eV between two con- 
secutive steps. Structural optimizations were performed 
using a conjugate gradient algorithm with a convergence 
criterion of 10 -4 eV/A. Pressures on the lattice unit cell 
are decreased to values less than 1.0 kBar. The adsorp- 
tion of a single CI atom to graphene surface is treated 
using the supercell geometry, where single CI is adsorbed 
to each (4x4) supercell. In the self-consistent potential 
and total energy calculations using a (4x4) supercell of 
chlorographene, a set of (13x13x1) k-point sampling is 
used for Brillouin zone (BZ) integration. The sampling 
of BZ is then scaled according to the sizes of the su- 
percells used for other systems. Ground state electronic 
structures are calculated by applying a dipole correction 
to eliminate the artificial electrostatic field between peri- 
odic supercells. For the charge transfer analysis, the ef- 
fective charge on atoms are obtained by Bader method^ 2 - 

The stabilities of structures having various CI cover- 
age are examined by the calculation of phonon frequen- 
cies for q-wave vectors over BZ by using both small dis- 
placement method (SDM)^ 3 . and density functional per- 
turbation theory (DFPT)4i We used 196-atom supercell 
of chlorographene, q-point sampling grid of 3x3x1 and 



0.01 A displacements in calculations using SDM. DFPT 
part of the phonon calculations were performed by us- 
ing 6x6x1 grid of q-points for chlorographene unitcell. 
A given structure is considered to be stable, if vibration 
frequencies are positive for all q-points in BZ. 

The energy band gap, which is usually underestimated 
by DFT, is corrected by frequency-dependent GWo 
calculations^ In GWo corrections screened Coulomb po- 
tential W, is kept fixed to initial DFT value Wo and 
Green's function G, is iterated four times. Finally, the 
band gap of CC1 is calculated by using (12x12x1) k-points 
in BZ, 20 A vacuum spacing, default cut-off potential for 
GWo, 160 bands and 64 grid points. 

The binding energy of a single adatom X to the 
graphene supercell (i.e. X=H hydrogenation; X=F fluori- 
nation; X=C1 chlorination) is calculated according to the 
expression, Eb = Er[Gr] + Et[X] — Ex[Gr + X], in terms 
of the ground state total energies 4 ^ of bare graphene 
Et [Gr] , free X atom Et [X] , and X adsorbed to graphene 
Et[Gt + X]. Accordingly, Eb > indicates a bound 
state. Similarly, the formation energy of a single adatom 
adsorbed to graphene supercell relative to X2 molecule 
is defined as E } = E T [X 2 }/2 + Et [Gr] - E T [Gr + X], 
where Et[X2\ is the ground state total energy 4 ^ of X2 
molecule. For an exothermic process Ej > 0. The 
cohesive energy of graphene fully covered with X, CX 
relative to free C and X atoms is defined as E co h = 
2E T [X] + 2E T [C] - E T [CX], where E T [CX] is the total 
ground state energy (per unit cell) of CX. By definition 
E co h > indicates binding relative to individual con- 
stituent atoms. 



III. ADSORPTION OF SINGLE CHLORINE 

Understanding of the adsorption process of a single CI 
adatom on graphene is essential for the investigation of 
its chlorinated derivatives. The Cl-Cl coupling between 
adsorbed CI atoms is crucial. Here, the CI coverage is 
defined as the ratio of the number of CI atoms Nci to 
the number of carbon atoms Nc in a supercell; namely 
G = Nci/Nq- A single CI atom adsorbed to graphene 
is represented by a system where a CI atom is adsorbed 
to each (4x4) supercell. This actually corresponds to 
a uniform coverage of 8=0.03125. Based on the anal- 
ysis discussed in the next section we found the size of 
(4x4) supercell is sufficient to neglect the Cl-Cl coupling. 
Even if the Pauling scale of electronegativity of specific 
atoms alone cannot provide criteria for the character and 
strength of the bonds between those atoms, it usually in- 
dicates useful correlations. For example, according to the 
Pauling scale of electronegativity (H: 2.20, C: 2.55, CI: 
3.16 and F: 3.98) the binding of CI atoms to graphene 
can be expected to be stronger than that of hydrogen 
adatoms (0.98 eV) and weaker than the fluorine adatoms 
(2.71 eV). Three different adsorption sites can be fore- 
seen for the adsorption of CI on graphene as described in 
Fig. [Ua) . These are the hollow (H) site above the cen- 
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FIG. 1: (Color online) (a) Variation of energy for a single CI 
adatom along the symmetry direction of a hexagon. Zero of 
energy is set to the energy of T-site. The diffusion path with 
the lowest energy barrier of Q=13 meV between two adjacent 
T-sites are marked with thick red/ dashed lines, (b) Energy 
landscape of a single CI adatom adsorbed to graphene. Dark 
(light) colors represents the top (hollow) sites, (c) Potential 
energy contour plots (paraboloid) of CI adatom around the T- 
site. The jump frequency of CI atom v, for different directions 
are calculated from this paraboloid, (d) Band structure of a 
single CI adsorbed to each (4x4) supercell of graphene and 
corresponding total (TDOS) and orbital decomposed (PDOS) 
densities of states. The zero of band energy is set to the Fermi 
level. 



ter of hexagon, the top (T) site on top of C atoms, the 
bridge (B) site above the middle of C-C bonds. Among 
these sites, the strongest binding of single CI to graphene 
is attained at T-site with a binding energy of £'6=1.10 eV 
(£ fc calculated using GGA+vdW 1.16 eV). Earlier stud- 
ies dealing with graphene-Cl interaction calculated the 



binding energies of CI adatom ranging 1.05 eV— using 
(4x4) supercell and 0.8 eV.— Using a methodology simi- 
lar to that used in the present work the binding energy of 
CI adatom is calculated to be 1.13 eV<2£ Applying correc- 
tion based on hybrid functional^ and hence performing 
GGA+HSE06 calculations we found the binding energy 
of CI as 0.7 eV. 

The variation of the total energy of an adsorbed CI 
atom at the symmetry points and along the symmetry 
(T-B-H-T) directions of a hexagon in the (4x4) supercell 
is shown in Fig.Q](a). At each point on the energy curve, 
x- and y-positions of adsorbed CI atom are fixed, its z- 
height, as well as positions of all C atoms in the (4x4) 
supercell are optimized by minimizing total energy and 
atomic forces. Using this energy curve we reveal the en- 
ergy barriers to be overcame by the adsorbate migrating 
or diffusing on the surface of graphene. The minimum 
energy barrier occurs at the B-site between two adjacent 
T-site. The relevant energy barriers are shown on the cal- 
culated energy landscaping presented in Fig. [T](b). This 
analysis suggests that a diffusing CI atom can take a path 
of minimum energy barrier following the edges of hexagon 
from T-site to T-site through the barrier of Q=13 meV 
at B-site. This energy barrier is in fair agreement cal- 
culated by Wehling et al^ We note that this barrier is 
very low and allows CI atom to migrate on the surface of 
graphene as if it rolls along the edges of hexagon. This 
is a remarkable situation, where the migration or diffu- 
sion of CI occurs with almost no barrier, but it remains 
bound to the surface. This property of CI adatom on 
graphene heralds a number of possible applications, such 
as superlow sliding friction, sensors and devices for en- 
ergy harvesting-^ The binding energy of an isolated CI2 
molecule is calculated to be 3.58 eV,— leading to negative 
formation energy Ef <0. Then the adsorption of chlo- 
rine on graphene is an endothermic process and hence it 
can be achieved by, for example, the assistance of light 
absorption.— 

In Fig. [T] (c) the energy plot Ej<(x,y) (or energy 
paraboloid) of CI adatom displacing from its equilibrium 
position at T-site is deduced approximately from the en- 
ergy landscaping shown in Fig. [T] (b). When displaced 
along T-H and T-B directions, the CI adatom undergoes 
different energy variations. Using this energy plot and 
the model related with simple harmonic oscillator we es- 
timate the values of jump frequencies of adatom towards 
H- and B-sites to be vt^b ~ 0.97 THz and vt^h ~ 
1.68 THz, respectively. The diffusion constant at a given 
temperature can be obtained from a simple expression,— 
D = va 2 exp(—Q / ksT) in terms of energy barrier Q, lat- 
tice constant a and the characteristic jump frequency v. 
Using these values the diffusion constant at T=300 K 
is calculated for the barrier Q=13 meV to be £>=3.57 
err? js. This indicates that CI adatom migrates readily 
over the graphene layer. 

The length of Cl-C bond is found to be d C /-c=2.54 A. 
This value is larger than the bond length calculated from 
the sum of covalent radii of CI and C atoms^ namely 
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Cl-Cl Distance (A) 

FIG. 2: (Color online) The interaction energy between two 
CI atoms adsorbed to the same side of a (6x6) supercell of 
graphene. The zero of energy is set to the energy of CI2 plus 
graphene. nn denotes the nearest neighbor in graphene lattice 

tci + r c = 1.75A, as well as the bond lengths of CCI4 
(1.77 A). This situation indicates that the character of 
this bond is rather different from C-Cl covalent bond in 
CCI4. When single CI is adsorbed to graphene, 0.44 
e is transferred from graphene to CI adatom according 
to Bader analysis^ The calculated charge transfer, as 
well as the analysis of charge density distribution indi- 
cates an ionic character for the Cl-C bond.— The neg- 
ative charge on CI induces a dipole moment of p=0.67 
eA when a single CI atom is adsorbed to each (4x4) su- 
percell of graphene.— Consequently, the work function 
of bare graphene which is 4.49 eV increases to 6.53 eV 
upon CI adsorption. This the energy required to remove 
an electron from the Fermi level to the vacuum energy 
at the side where CI atoms are adsorbed. Incidentally, 
while the atomic configuration of graphene beneath CI 
adatom is rather flat, H and F adatoms on graphene im- 
pose buckling resulting in a covalent character i 14 ' 23 As 
we will discuss in the next section, charge transfer from 
C to CI and hence the character of the bonding in two 
sided adsorption is dramatically different. 

While single CI adatom on the surface of graphene can 
be viewed as an impurity leading to localized (or reso- 
nance) states, our model representing the adsorption us- 
ing a (4x4) supercell gives rise to an energy band struc- 
ture, where dispersion of bands related with adsorbate 
can be taken as the measure of adsorbate-adsorbate cou- 
pling. The band structure of a single CI adsorbed to each 
(4x4) supercell of graphene presented in Fig. [T] (d) has 
filled chlorine states appearing as flat bands. In view of 
the negligible dispersions of these bands associated with 
CI and the analysis of the bands corresponding to (5x5) 
and (6x6) supercells in the following section we concluded 
that a CI adsorbed to each (4x4) supercell can mimic suc- 
cessfully the single, i. e. isolated CI adsorbed to graphene 
surface. Because of 7 valence electrons and hence a sin- 



gle unfilled 3p orbital occurring below the Fermi level Ep 
of bare graphene, CI adsorbed graphene is metalized and 
becomes magnetic with /x=0.56 ijlb due to broken spin de- 
generacy. The linearly crossing bands of bare graphene 
is raised above Ep. The orbital decomposed density of 
states in Fig. Q] (e) show that the flat bands below Ep 
are derived mainly from Cl-3p x and C\-3p y orbitals. 



IV. COVERAGE OF GRAPHENE BY 
CHLORINE ADATOMS 

The interaction between two adsorbed CI BjS ct func- 
tion of their separation is important to understand the 
coverage dependent properties and stability of CI covered 
graphene. To this end we consider two CI atoms adsorbed 
to the (6x6) supercell of graphene and calculate the total 
energies as a function of the separation between them. 
Chlorine atoms are placed at specific adsorption sites by 
fixing their x- and y-coordinates, but their z-coordinates, 
as well as the positions of other C atoms in the supercell 
are relaxed. In Fig. [2] we show how the Cl-Cl interac- 
tion changes with their separation. The Cl-Cl coupling 
is practically negligible when their separation greater or 
equal to 3.6 A corresponding to fourth nearest neighbor. 
For smaller separations corresponding to third nearest 
neighbor separation of graphene with the threshold dis- 
tance, they form CI2 molecule and desorb from graphene 
surface, since Cl-Cl interaction energy becomes stronger 
than the sum of the binding energies of two CI atoms 
to graphene (2.20 eV). Accordingly the gain of energy 
through the formation of CI2 is 1.26 eV per molecule. 
This explains also why chlorination of graphene is en- 
dothermic. There is a weak vdW interaction between 
bare graphene and CI2 molecule. The maximum interac- 
tion of 144 meV occurs when the molecule is parallel to 
a C-C bond and 3 A above it. 

Next we examine how the electronic structure changes 
with one sided uniform coverage. In Fig. [3] we show the 
band structures calculated for single CI atom adsorbed 
to (mm) supercell for n=2,3,5 and 6. Here n=2 and n=3 
correspond to 0=0.125 and 0.056. The phonon calcula- 
tions for these two coverage values are found stable, since 
the frequencies of phonon modes are positive for any k- 
point in BZ. The smallest separation between CI atoms is 
larger than the threshold distance of 2.8 A (i.e. the third 
nearest neighbor distance) even for 9=0.125. We see that 
with increasing 9 the dispersion of CI bands increases 
slightly and the linearly crossing bands start to split. The 
magnetic moment also increases with coverage, since the 
splitting of spin- up and spin-down bands increases. We 
found /j,=0.64 fj,B and //=0.56 fiB for 9=0.0555 (or n=3) 
and 9=0.020 (or n=5), respectively. On the other hand, 
a single CI adsorption to each (lxl) unit cell of graphene 
(which corresponds to 9=0.5) is unstable due to Cl-Cl 
distance smaller than 2.8 A. Since the uniformly chlori- 
nated graphene is found to be metallic for any coverage 
studied here, the experimental studies measuring a band 



5 




FIG. 3: (Color online) Energy band structure of a single 
CI atom adsorbed to each (nxn) supercell of graphene for 
n=2,3,5 and 6, which correspond to the one-sided uniform 
coverage G = l/2n 2 . For n > 2 the distance between nearest 
CI adatoms is larger than the the threshold distance described 
in Fig. [2] Whereas stable CI coverage with n = 1 (or ©=0.5) 
cannot be achieved due to the strong Cl-Cl coupling. The 
units of magnetic moments p is Bohr magneton per (nxn) 
supercell. The zero of energy is set at the Fermi level, Ef- 

gap of ^0.05 eV should correspond to two sided coverage. 

To understand the formation of chlorinated domains 
on graphene, and hence to analyze the kinetics of CI cov- 
erage, we next examine the adsorption of second CI atom 
at different sites at the close proximity of the first one in 
(4x4) supercell. When a CI is placed on top of a carbon 
atom, second CI can be adsorbed on six relevant sites on 
the hexagonal carbon ring as shown in Fig. 2] (a) . 

Among the various possible cases, the ortho top- 
bottom arrangement is the most favorable one. While 
adsorption of a single CI atom yields 0.05 A buckling of 
graphene lattice, adjacent C atoms forming C-Cl bonds 
in ortho top-bottom configuration are buckled by 0.46 A. 
In Fig. HJb) the contour plots of total charge density on 
the planes passing through the bonds clearly explain the 
ionic character of C-Cl bond of a single CI adatom. In the 
case of ortho top-bottom the bonding between CI adatom 
and C atom of graphene has sp 3 -type covalent character. 
While the charge density between CI and C atoms is very 
low in the ionic case, the bond charge is enhanced in the 
covalent C-Cl bond. The binding energy of covalent bond 
is stronger (£^=1.53 eV) and consequently its length is 
shorter (rf=1.90 A, which is comparable with the empir- 
ical covalent radii of the C-Cl bond or the bond length 
dc-ci=l-77 A of CCI4 molecule). The binding ener- 
gies per bond in the case of meta and para top-bottom 
configurations are i?b=1.10 eV and £^=1.31 eV, respec- 
tively. In spite of the fact that the ortho top-bottom 
configuration is 0.52 eV less energetic than the forma- 
tion of CI2 molecule, it can be stable since the formation 
of CI2 molecule from two CI atoms at different sides are 
hindered. In para and meta top-bottom structures, the 
binding energies of two bonds are further lowered from 
the formation energy of CI2 molecule. 

The situation is however dramatically different for one 
sided (top-top) configurations presented in the second 
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FIG. 4: (Color online) (a) The atomic structure two CI atoms 
adsorbed to a (4x4) supercell of graphene. In three differ- 
ent configuration illustrated by top panels, namely ortho top- 
bottom, para top-bottom and meta top-bottom, two adsorbed 
CI atoms are stable. AE indicates their energies relative 
to the total energy of the ortho top-bottom configuration. 
Double sided adsorption imposes a local buckling in planar 
graphene. Three one-sided configurations, ortho top-top, para 
top- top and meta top-top are not allowed; CI atoms cannot 
be bound to graphene, they rather form CI2 molecule. Large 
green and small gay balls represent CI and C atoms, respec- 
tively, (b) Contour plots of the total charge density of a single 
Cl-C bond and two Cl-C bonds in ortho top-bottom config- 
uration. Contours spacings between 0.025e/A 3 and 1.0 e/A 3 
are 0.025 e/A 3 . 



row of Fig. [4] (a). For these three one-sided configuration, 
namely ortho, para and meta top-top configurations the 
Cl-Cl distances are within the third nearest neighbors. 
Owing to the strong Cl-Cl coupling, the formation of 
CI2 molecule is energetically favored. Accordingly, unlike 
the half fluorinated graphene^ (C2F), one-sided densely 
chlorinated graphene with the coverage 0=0.5, namely 
C2CI, cannot be realized. 

We now concentrate on the chemical conversion of 
graphene to a fully chlorinated graphene structure 
(0=1.0), called chlorographene. In addition to well- 
known boat, chair and nonbonding chair conformations 
one can also consider zigzag and armchair stoichiomet- 
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(a) Boat Chair Zigzag Armchair 




unstable stable unstable unstable 




FIG. 5: (Color online) (a) Top view of atomic structures of 
boat, chair (i.e. covalently bonded and buckled graphene 
layer is sandwiched between two planar CI layers) , zigzag and 
armchair conformations. Green and small gay balls represent 
CI and C atoms, respectively, (b) Side view of nonbonding 
chair conformation consisting of one planar graphene layer 
sandwiched by two planar CI layers and its calculated phonon 
dispersion curves. Low frequency phonon modes shown by red 
lines are related to adsorbed CI atoms. These modes have 
imaginary frequencies and hence they are unstable. 



ric chlorographene configurations shown in Fig. [5] Here, 
boat, chair, nonbonding chair and zigzag conformations 
are treated using (2x2) supercell; a (4x4) supercell is re- 
quired for the armchair conformations. As a consequence 
of the strong Cl-Cl interaction, boat, zigzag and arm- 
chair configurations of CI atoms cannot remain stable on 
graphene. 

Here the nonbonding chair conformations in Fig. [5] (b) 
consists of a planar graphene layer sandwiched between 
two CI layers deserves a detailed discussion. In this con- 
formation CI atoms are placed to the alternating T-sites 
at both sides of graphene. Earlier theoretical studie a 30 i 31 
have predicted that this metallic structure has highest co- 
hesive energy among conformations described in Fig. [5] 
Present calculations predict £/r=23.35 eV/per primitive 
cell and cohesive energy E co h=2.79 eV / per primitive cell. 
Even if the structure optimization using conjugate gra- 
dient method performed in the (lxl) hexagonal unit cell 
finds the nonbonding chair structure stable, we carried 
out an extensive analysis of this conformation. We found 
that the attractive interaction between each CI layer and 
underlying graphene is rather small (50 meV per unit 
cell) in spite of the fact that graphene lattice expanded 
by 2.8% and the system attained a total energy <~3.31 
eV higher than that of graphene. The weak interaction 
between CI and C atoms explains why the distance from 



CI layer to graphene is large (3.51 A). This large distance 
and minute charge transfer from carbon to CI atom can 
be contrasted with the C-Cl distance (2.54 A) and the 
transfer of 0.44 electrons from C to CI of the ionic bond 
of single CI atom adsorbed to graphene as discussed in 
Sec. III. This paradoxical situation occurred due to the 
strong intra-layer coupling among CI atoms in both Cl- 
layers at both sides of graphene. In fact, the binding 
energy of CI atoms of a single CI layer having the same 
atomic configuration as the CI layers of nonbonding chair 
conformation is found to be 1.41 eV per CI atom. It 
appears that high cohesive energy of nonbonding chair 
conformation relative to free graphene and free CI atoms 
is attained mainly by the cohesion between CI atoms, 
but not by the interaction between CI and carbon atoms. 
In addition, GGA+vdW calculations performed for non- 
bonding conformation resulted with optimized structure 
close to that obtained by using LDA. 

Under these paradoxical situations we further exam- 
ined the stability of the metallic nonbonding chair con- 
former by carrying out ab-initio calculations of phonon 
frequencies in the BZ. As seen in Fig. [5] (b), the cal- 
culated phonon branches, which are associated with CI 
in-plane and out-of-plane modes and hence are practi- 
cally isolated from graphene modes, have imaginary fre- 
quencies near M- and A'-high symmetry points. This 
clearly indicates that nonbonding chair structure is un- 
stable at T=0 K. In fact, a minute displacements or per- 
turbation of atoms in CI plane is resulted in the break- 
down of Cl-layers, if nonbonding chair structure is treated 
using a (4x4) supercell. It appears that structure opti- 
mization in earlier studie o 30 ' 31 carried out in a single cell 
has limited the degree of freedom of CI atoms preventing 
them from reconstructions involving large displacements 
or from forming molecules. Accordingly, the stability im- 
posed by the structure optimization using primitive cell 
was unrealistic. Even though the nonbonding chair con- 
formation has the lowest energy among other conformers 
including (covalently bonded) chair structure in Fig. [SJ it 
corresponds either to a very shallow local minimum or to 
a saddle point for specific directions of atomic displace- 
ments in Born-Oppenheimer surface. 



V. STABLE FULLY CHLORINATED 
GRAPHENE: CHLOROGRAPHENE 

A. Structural Properties 

In contrast to nonbonding chair conformation, our 
analysis indicates that like graphane (CH) and fiuo- 
rographene (CF), the chair structure, where one chlo- 
rine atom is attached to each carbon atom of buckled 
graphene alternatingly from top and bottom sides can 
be a stable structure. The optimized atomic structure 
and the hexagonal lattice are shown in Fig. |H] (a). The 
present study predicts its total energy £r=22.24 eV/per 
primitive cell and its cohesive energy £' co ^=1.68 eV. In 
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TABLE I: Calculated values for graphene, and graphene derivatives, such as graphane CH, fluorographene CF and chloro- 
graphene CC1. These are lattice constant (a); C-C bond distance (dec); C-X bond distance (dcx)', thickness of the layer (t); 
photoelectric threshold (<£>); charge transfer from C to X (Ap); cohesive energy per unit cell (E co h), formation energy (Ef), 
desorption energy of a single X atom from the CX surface (Ed ea ); direct band gap (E g ); band gap corrected with GW D , _E^ W ° 
; in-plane stiffness (C); Raman active modes. Energies are calculated in (2x2) supercell. 
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the rest of paper we name this stable structure as chloro- 
graphene or CC1, which should be distinguished from the 
unstable, nonbonding chair structure presented in the 
previous section in Fig. [5] (b). The remarkable situa- 
tion is that the sp 2 -bonded planar graphene is buckled by 
(5=0.50 A, once it is converted to CC1. This way the local 
configuration of three folded sp 2 -bonding has changed to 
four folded sp 3 -bonding, reminiscent of tetrahedrally co- 
ordinated diamond structures. Formation of four folded 
sp 3 -bonding through buckling is essential for the stabi- 
lization of CC1, despite the C-C distance increased from 
1.42 A to 1.72 A. On the other hand, while the C-Cl ionic 
bond in the adsorption of single CI adatom is 2.54 A in 
Sec. Ill, it transforms to a covalent bond by contracting to 
1.73 A in chlorographene, which is also very close to the 
bond length of CCU molecule. The bond length of chloro- 
graphene can also be contrasted with the Cl-graphene 
distance in nonbonding chair structure discussed in the 
previous section. We note that chlorographene is stable 
despite its cohesive energy is smaller than that of the un- 
stable nonbonding chair conformation. This is due the 
fact that chlorographene corresponds to a local minimum 
in Born-Oppenheimer surface. 

In Table U we compare the structural parameters, rele- 
vant energies like cohesive energy, formation energy and 
desorption energy of graphene, CH, CF and CC1, all cal- 
culated by using LDA. It is seen that chlorographene, 
CC1 has binding energy E u =l.68 eV per unit cell rela- 
tive to bare graphene and free CI atom, and lowest E co h 
(19.60 eV) per unit cell among other possible graphene 
derivatives. However, once CC1 is synthesized a desorp- 
tion energy of Edes =1-28 eV is required to take a single 
CI atom out. 



B. Vibrational Properties and Raman Spectra 

Chlorographene (CC1) posses D 3( z point group symme- 
try. Phonon bands of CC1 is shown in Fig. [6jb). It is 
seen that all the phonon modes have positive frequen- 
cies and hence the predicted structure of chlorographene 
is stable at T=0 K. Unlike the phonon bands of non- 
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FIG. 6: (Color online) (a) Top, side and tilted views for the 
atomic structure of chlorographene CC1 layer having hexag- 
onal lattice and honeycomb structure. Carbon and chlo- 
rine atoms are indicated by gray (dark) and green (light) 
balls, respectively. Calculated structural parameters are indi- 
cated, (b) Phonon bands of chlorographene calculated using 
by SDMi2 and DFPT— methods. 



bonding conformer, C-Cl bonds of chlorographene par- 
ticipate to all acoustic phonon branches. Therefore, the 
phonon spectrum of chlorographene differs from the spec- 
trum of graphene. Though the chlorographene belongs to 
the same space group with CH and CF ) 22 i 23 the phonon 
frequencies are lowered (softened) due to the saturation 
of C atoms with heavy CI atoms. While LA and TA 
modes have linearly dependent to phonon wave vector, 
ZA mode has a quadratic dispersion in the vicinity of P- 
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FIG. 7: (Color online) (a) Electronic band structure of chloro- 
graphene CC1. The band gap is shaded yellow. The GW„ cor- 
rected valance and conduction bands are shown by dashed line 
and red balls. The zero of energy set to the Fermi level Ep. 
(b) Density of states projected to various orbitals (PDOS). 



pointi 22 ' 23 i 52 As compared to single layer graphene, LO 
and TO (ZO) optical modes are softened from 1600 (900) 
to 1061 (421) cm -1 due the existence of surrounding CI 
layers. In general, both phonon bands calculated using 
SDM and DFPT agree well, but they slightly differ at 
M and K point for optical phonon branches near at 500 
cm" 1 . 

Group theory analysis shows that the decomposition 
of the vibration representation at the T-point is T = 
2A\ g + 2A 2u + ^E g + AE U . Among these, the modes at 
105, 398, 715 and 1042 cm" 1 are bond stretching modes 
and are Raman-active. Raman mode A\ g at 1042 (398) 
is entirely due to the out-of-plane vibration of C and CI 
atoms moving in the same (opposite) direction with re- 
spect to each other. The observation of these Raman ac- 
tive modes are expected to shed light on the CI coverage 
and the structure of chlorinated graphene. The obser- 
vation of characteristic D-peak at 1330, G-peak at 1587 
and 2D-peak at 2654 cm" 1 from chlorinated graphene in- 
dicates low coverage of Cl^ Since the Raman measure- 
ments were performed in the range of 1250-3500 cm" 1 , 
possible Raman-active peaks originating from chlorine 
atoms could not be observed. 



C. Stability of CC1 at finite temperature 

Here we investigate whether CC1 is stable at finite tem- 
perature, even though all phonon modes having positive 
frequencies in the Brillouin zone indicate its stability at 
T=0 K. This is achieved by carrying out ab-initio molec- 
ular dynamics (MD) calculations at finite temperatures 
using (6x6) supercell involving 144 C and one CI atoms. 
Here we summarize our findings. At relatively low tem- 
peratures, for example at T=500 K, the perfect CC1 re- 
mained stable even after ^6000 time steps of 2 fs. Even 
if 6000 time steps large for ab-initio calculations, but low 



to attain a reliable statistics, this result suggests the sta- 
bility of perfect CC1 near room temperature. However, 
CI atoms dissociate from CC1 at 1000 K. This may be 
interpreted that CC1 cannot be stable at elevated tem- 
peratures. Ab-initio MD calculations with single CI va- 
cancy at one side brought about the formation of a sec- 
ond vacancy at the other side after 200 time steps at 
500 K. Thereafter, the structure continued to be stable. 
However, an island of CI at both sides of graphene, for 
example, 6 CI atoms adsorbed to a hexagon of graphene 
alternatingly from the top and bottom side of graphene 
is found unstable. 

Our analysis using finite temperature ab-initio MD cal- 
culations led us to draw the following conclusion from 
the above results: A perfect CC1 is stable at T=0 K and 
can remain stable possibly at room temperature. The 
creation of a single vacancy at one site imposes the for- 
mation of a second vacancy at the opposite side. This 
pair of vacancy can survive at room temperature. How- 
ever, CC1 having vacancies or holes is vulnerable to dis- 
sociation through the formation of CI2 molecules. The 
negative formation energy underlies these instabilities. 



D. Electronic Properties 

We present the band structure and the orbital de- 
composed density of states of fully chlorinated graphene 
CC1 in Fig. [7] As a consequence of transition from 
sp 2 - to of sp 3 -type bonding and the mixing with CI 
orbitals through Cl-C covalent bonds, linearly crossing 
semimetallic bands of graphene changes dramatically and 
turn to a nonmagnetic semiconductor with a direct band 
gap of E g =1.21 eV (GGA+vdW value: 1.55 eV) at the 
center of BZ. The two-fold degenerate bands at the top 
of valence band are mainly composed of the p xy valence 
orbitals of CI and C atoms and belong to the E g irre- 
ducible representation. However, non-degenerate con- 
duction band edge is formed by the hybridization of C- 
Vxy-, C-pz, C-s and Cl-p z and belongs to the A 2u irre- 
ducible representation. Calculated values relevant for the 
electronic properties of CC1 are given together with those 
of graphene, CH and CF in Table Q] CC1 has the low- 
est band gap (1.21 eV) among graphene derivatives, i.e. 
graphane (E g =3.42 eV) and fluorographene (E g =2.96 
eV) . Our calculations also reveals that the experimentally 
observed energy gap of 0.045 eV may arise from the par- 
tially chlorinated regions. Since DFT usually underesti- 
mates the band gaps, the band gap of 3D bulk crystals are 
successfully corrected by GW Q self-energy method. In the 
present paper we apply GW method to correct valance 
and conduction bands and hence the band gap of CC1. 
We found the corrected band gap 4.33 eV. This is a dra- 
matic increase. Similar situation occurred for the correc- 
tion of band gaps carried out for single layer structures, 
such as BN2 3 ., CF2 3 . and MoS^. In particular, while 
LDA predicted band gap of single layer M0S2 agrees with 
experiment, GW Q correction yielded very large band gap. 
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This situation led us to question whether GW correction 
is suitable for single layer structures. 

While graphane have a positively charged surface due 
to the electron transferred from H to C atoms, CI lay- 
ers of CC1 are negatively charged since 0.13 electrons 
(GGA+vdW value: 0.10 electrons) are transferred from 
C to CI atoms. We also note that the effective charge 
is reduced from 0.42 to 0.13 electrons by going from 
ionic C-Cl bond in single CI adsorption to the cova- 
lent C-Cl bond of chlorographene. Since CC1 surfaces 
are negatively charged, it is possible to lower the photo- 
electric threshold of graphene from 4.49 eV to 3.67 eV 
by covering its surfaces with chlorine atoms. In con- 
trast, as one notices in Table U both fluorinated^ 3 - and 
hydrogenate d 14 ! 22 derivatives of graphene have photo- 
electric thresholds higher than graphene and CC1. The 
lower photoelectric threshold provide materials highly 
emissive surfaces. This facile photoemission feature is 
desirable for fast laser applications 



E. Mechanical Properties 

Earlier studies have shown that graphan o 15 i 16 and 
fluorographen o 22 ! 23 derivatives are strong materials like 
graphene. Among the three graphene derivatives chloro- 
graphene is the thickest one and we can expect some 
different mechanical properties. We will discuss the 
structural rigidity of chlorographene within the harmonic 
range of the elastic deformation, where the structure re- 
sponded to strain e linearly and reversibly. In the elastic 
range, in-plane stiffness C would be a good measure of 
the response of material. Here we use the expression 
to calculate in-plane stiffness, C = (1/Aq) ■ (d 2 Es/de 2 ), 
where Aq is the equilibrium area of the supercell and the 
strain energy is defined as the total energy at a given uni- 
form strain minus the total energy at zero strain, namely 
E s = Et{() — Et{c = 0). Here the strain in one direction 
is e = Ac/c i where cq is the equilibrium lattice constant 
of supercell and Ac is its stretching. For CC1 we calcu- 
lated C as 186 J/m 2 which identifies it as strong as CH 
and CF. 

In Fig. Ufa) we show the variation of the strain energy 
E s and its derivative with respect to the applied uniform 
strain, dE s /de. The latter is linear for small e in the har- 
monic range. The elastic deformation occurs until the 
maximum of dE s /de, whereby the structure attains its 
initial state when the applied strain is lifted. Beyond the 
maximum the structural instabilities sets in with irre- 
versible deformations. The region beyond the maximum 
is called plastic region. It seen that CC1 undergoes insta- 
bility, under 0.15 expansion of the lattice. Whereas, the 
critical expansion values are 0.20 and 0.24 for CF and 
CH, respectively. 

We also calculated the effect of elastic strain on the 
band gap of CC1, CF and CH and presented our results 
in Fig. [5{b). CH and CC1 have similar response to elastic 
strain: The increase of their band gaps for small strain is 




0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 



Ac/c„ Ac/c„ 

FIG. 8: (Color online) (a) The variation of the strain energy 
E s (curves at the right hand side) and its derivative with 
respect to applied uniform strain e, i.e. dE a /de, (curves at 
the left hand side) are calculated for CC1, CH and CF. After 
the maxima these structures become unstable and undergoes 
a plastic deformation, (b) Variation of the band gap with 
uniform strain. Calculations performed in (5x5) supercells. 



followed by a rapid decrease for large strain. In contrast, 
initially the band gap of CF does not show any signifi- 
cant increase with increasing strain; for small strain it is 
almost unaltered, but decreases rapidly for large strain. 
Because of its smaller band gap, the metalization of CC1 
occurs earlier in the course of expansion. 



VI. CONCLUSIONS 

Motivated with the recent work by Li et ali 2 ^, 
who achieved the photochemical chlorination of a bare 
graphene, we performed a first principle study on the 
chlorination of graphene starting from a single CI adatom 
adsorption to fully chlorinated graphene CC1. We found 
that even if a CI atom can be bound to graphene with a 
significant binding energy, it can migrate on the surface 
of graphene almost without an energy barrier. Formation 
of a CI2 molecule from two individual, migrating adatoms 
is energetically favorable when they are at close proxim- 
ity. Therefore the formation energy of the adsorption of a 
single CI atom relative to CI2 molecule is negative. In this 
respect, all possible configurations or decoration obtained 
by the chlorination of one surface of graphene and con- 
sidered in this study are not stable. On the other hand, 
the configuration, where two CI atoms are adsorbed from 
opposite sides of graphene to two adjacent carbon atoms 
forming sp 3 -type covalent C-Cl bonds and hence induc- 
ing buckling are stable, since the dissociation of these two 
CI atoms by forming CI2 molecule is hindered. Once the 
bare graphene is fully chlorinated from both side induc- 
ing buckling of carbon atoms, the resulting conformation 
named as chlorographene, is stable at T=0 K and also 
possibly at room temperature. With its 1.21 eV direct 
band gap, stiff mechanical properties and response to ho- 
mogeneous strain this material displays interesting prop- 
erties for future technological applications. On the other 
hand, the nonbonding chair conformation has cohesive 
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energy higher than that of chlorographene in the chair 
conformation and is found to be unstable. While both 
structures are not in ground state, chlorographene cor- 
responds to a local minimum on the Born-Oppenheimer 
surface. Of course, a state having the energy lower than 
those of both conformation is CI2 molecules, which are 
weakly attached to bare graphene. 
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